Theory of measurements of electrodynamic properties in anisotropic superconductors 
in tilted magnetic fields. Part II: high frequency regimes. 
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The model for high frequency electrodynamics in anisotropic type-II superconductors in the vortex 
state is studied considering arbitrary orientations between the applied field, the applied current and 
the anisotropy axis. An anisotropic treatment is provided for the vortex dynamics, taking into 
account all the phenomena relevant at high frequency, which include flux flow, pinning and creep. 
The coupling between vortex motion and high frequency currents is included, providing an entirely 
tensor model of the electromagnetic response to high frequency fields. Examples of data analysis of 
angular measurements are presented, showing how to derive the angular dependence of the material 
properties from the measured anisotropic response. Finally, the expression of the measured angle- 
dependent surface impedance in the largely used thin film geometry is computed. 
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I. INTRODUCTION 

Material anisotropy, which characterizes many superconductors of wide interest, 1 has a profound impact, among 
the others, on their vortex dynamics and on the related pinning phenomena. 

In this paper I focus on the high frequency regime of the electrical transport properties in the mixed state, largely 
studied because of the great deal of information that they can provide. 

Due to the anisotropy of the superconductor, the transport properties depend on the various angles between 
the anisotropy axes and the applied field and current, so that the measured quantities have a non-straightforward 
relationship with the material properties. 2 

While some aspects of the problem in the d.c. and (low frequency) a.c. regimes have been addressed in previous 
works, 3-6 I proposed 2 a generalized treatment centered on the force equation for the vortex motion in the linear regime, 
including both material anisotropy and pinning, in uniaxial anisotropic superconductors in a magnetic field applied 
with generic orientation. The very different free flux flow regime (dominated by dissipation) and pinned Campbell 
regime (dominated by pinning) were addressed. The tensor expressions for the resistivity and for vortex parameters 
like the viscous drag, the vortex mobility and the pinning constant were given, and the measurable quantities for 
arbitrary angles between magnetic field, current and the anisotropy axis were derived. 

In this work I extend the treatment of Rcf. 2 to the high frequency regime, where additional phenomena emerge. 
Indeed, by increasing the frequency, dissipation and pinning effects become comparable, requiring to be simultaneously 
taken into account. Moreover, even within the limit of small currents, vortex creep effects become relevant. 7-11 Finally, 
the coupling between vortex motion and the high frequency currents (due to both superfluid and quasiparticles) 
contribute significantly. 

This work is organized as follows: first, I generalize the anisotropic vortex dynamics model to include both dissipa- 
tion and pinning effects, the latter including the effects of thermal depinning/creep (Sec. II). Contextually I provide 
examples of experimental data analysis (Sec. III). Second, I consider the full coupling between vortex motion and 
the high frequency currents (Sec. IV A). Finally, I provide an application example of the full model by computing 
the expression of the measured surface impedance in the largely used thin film geometry (Sec. IV C). 



II. A.C. VORTEX MOTION RESISTIVITY 

In this Section I address the issue of the vortex motion resistivity tensor p v in a uniaxial superconductor, and 
I calculate explicit expressions for the angular dependence, with arbitrary angle between the applied field B, the 
applied current density J and the crystal axis. In order to introduce all the relevant quantities, I first recall the 
general models describing the high frequency vortex dynamics with reference to the commonly studied configuration 
of isotropic superconductors with J _L B (Sec. II A). Second, I develop the full treatment for uniaxially anisotropic 
superconductors (Sec. II B). Finally, some examples of data analysis based on the obtained results will be illustrated 
(Sec. III). 



A. Short review of scalar models 

The scalar force equation describing the vortex motion in an isotropic superconductor with isotropic (point) pinning 
and J _L B is, in the sinusoidal regime e lwt : 12 ' 13 

k 

T]V+ t^-V = $ J + F therm (1) 

where v is the vortex velocity, rj is the vortex viscosity, k p is the pinning constant (also called the Labusch parameter) 
yielding an elastic recall force —k p u, with u = v/(\w) the vortex displacement from the pinning center, and F t herm is 
a stochastic thermal force causing thermal depinning (vortex creep). Different approaches, 7 ' 8 with different ranges of 
applicability, 14 have been proposed to take into account creep effects. As an illustration, I follow here the description 
of thermal depinning in terms of the relaxation of the pinning constant k p (t) — k p e~ l l TtH (Rcf. 7). The characteristic 
time for thermal activated depinning is: 



(2) 
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where t p is the inverse of the (dc)pinning angular frequency r p 1 = lo p = rj/k p (which will be commented on later) 
and U is the activation energy. Equation (1) can then be rewritten as: 7 

k 1 

TjV+T?-- —V = T] C V = $ J (3) 

ID 1 — 

where the complex viscosity rj c has been introduced: 

*-*('-^T^r:) W 

The corresponding scalar vortex motion resistivity is: 

_<M? _$ B £ , + i^ 

^~~~ITT^ (5) 



where the characteristic angular frequency wo is: 
and the creep factor e' is: 



W = T^ 1 + T p 1 (6) 



e' = — (7) 

1 + c k bt 



For {/ — > oo the creep is negligible, e' — > and u; — > w p ; consequently the vortex motion resistivity becomes: 

P*= ~ rr^E = [Pff - l Pc ) ( g ) 

where pg = $oB/rj is the flux flow resistivity and pc = $o-B/(fc p w) is the well-known Campbell resistivity. 2,15 

This limit corresponds to the Gittleman-Rosenblum (GR) model. 12 From Eq. (8) it can be seen that the pinning 
angular frequency uj p marks the transition between a "low frequency" and a "high frequency" regime: for w C w P 
the pinning force dominates over the viscous drag, yielding p v — > ipc, while for lo >• cj p , a purely dissipative flux flow 
regime is recovered with p v « pg yielding the same behaviour as in d.c. with no pinning. 

Before concluding this short review, it is worth recalling the definition of the often used 16-18 dimensionless ratio r: 

- - ^ (9) 



which, if creep is negligible (GR limit) yields: 



R(p„) 



r = ^ = Pff (10) 

LO p C 



The r parameter can be directly computed from the complex resistivity p v and it is unaffected by any systematic error 
in the experiments. Physically, it allows to easily evaluate whether the vortex dynamics is in the pinning dominated 
(r ^> 1) or flux flow dominated (r <C 1) regime. 



B. The anisotropic model 

In this Subsection I introduce the anisotropy following the approach previously used. 19 All symbols will follow the 
notation introduced in Part I. 2 For the sake of completeness, I recall here some basic features. The crystallographic 
axes arc taken as x = a, y = b and z = c. The latter is the axis of anisotropy. Figure 1 shows the chosen frame of 
reference and the magnetic field B = BB = B(s'm8 sm<j), sin 9 cos 4>, cosO) = (B x , B y , B z ) with generic orientation. 
The anisotropy is considered through the mass tensor: 19 



(m ab \ (l \ 

m ab = m ab 1 = mM 
\ m c ) \0 7 2 y 



(11) 
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FIG. 1. Principal frame of reference. The magnetic induction field B is also depicted, applied along a generic direction at the 
polar 6 and azimuthal (j> angles. 



having denned the in-plane mass m a \, — m, the out-of- plane mass m c and the anisotropy factor j 2 — m c /m a b. The 
notations A, A and A denote a tensor/matrix, a vector and a unit vector, respectively. A static magnetic field is 
uniformly applied to the superconductor with generic orientation. Vortices are considered as straight and rigid lines, 
so that more complex phenomena such as helical instabilities of the vortex lines 22 are not taken into account. I assume 
the validity of the London approximation B ~ p H. 

The force equation in anisotropic superconductors including both viscous drag and pinning, momentarily neglecting 
vortex creep (which will be treated later), can be written as: 

fjv + ^-k p v = $ Jx B (12) 

1LU 



which, by defining the complex viscosity tensor as: 



(13) 



can be also written as: 



fj c v = $ JxB (14) 



It is evident that the force equation (14) is formally equivalent to force equations written for the flux flow and 
Campbell regimes in Ref. 2. Hence, it is straightforward to apply the electrodynamics model developed there 
to the scries of quantities [fj c , fjc^ , p~v\ crtf\ p, v , p v ] which are equivalent to the series [fj, ff % \ Pg\ &ff\ fiv, Pff] an d 

[kp 1 ^ /(iuj),ip^) ,ia-Q\fl v ,ipc}, previously extensively studied. 2 It should be recalled that an essential feature of the 
anisotropic electrodynamics resides in the difference between measurable properties and material properties. As an 
example, the measured resistivity tensor p is different from the intrinsic resistivity tensor . Here and henceforth, 
the superscript "(i)" indicates material, intrinsic properties. 
As a first result, within this framework one can write: 2 

fh = -B x f$ ) B x (15) 

where (see Appendix of Ref. 2): 

S v := 




Equation (15) and (13) allow to write the "intrinsic" tensor: 



#=^)-i^ (is) 



In Ref. 2 I obtained, neglecting Hall contributions and for weak random point pinning: 

fj^(B,e)= V ^(B/B c2 (0))M- 1 (17a) 
%P(B,6) - k^B/B^M- 1 (17b) 
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where the relation rj^ = $0-6/ ' pf} connects the viscous drag and the flux flow resistivity intrinsic tensors. From 
Eqs. (17), it is seen that 77W and k p ^ obey the so-called angular scaling law 4,23 ' 24 , according to which, in the London 
approximation, a thermodynamic or intrinsic transport property of a uniaxially anisotropic superconductor depends 
on the applied field magnitude B and directional angle 9 only through the ratio B/B c2 {9). Substituting Eqs. (17) 
into (16) yields: 




> P {B/B c2 {9)) 



fjW(B/B c2 (9)) (18) 

It is evident that fj^ inherits the angular dependencies and anisotropic properties from fj^ and k p \ represented by 
M and by the scaling law. 

I stress that, in obtaining Eq. (18), I have used a very important result for the analysis of the experiment, namely: 



k p %{B/B c2 (9)) 
V ^(B/B c2 {9)) 



= uj Ptii (B/B c2 (9))=u p (B/B c2 (9)) (19) 



which holds for i = 1..3, i.e. for all the principal axes directions. The above equation shows that, contrary to the 
viscosity or the pinning constant, the pinning frequency is a scalar. 25 This fact implies that, whichever orientation is 
set for B and J, the vortex system will always have the same pinning frequency at fixed B / B c2 {9). This property 
suggest a straightforward method to check whether directional defects influence vortex motion: by plotting uj p vs the 
scaled field B/B c2 (9), any deviation from a scaling curve should indicate the influence of some directional effect, other 
than the material anisotropy. Obviously, if extended defects are present as sources of pinning, Eq. (17b) does not 
hold and the above no longer applies. 

The intrinsic conductivity and resistivity tensors are: 2 



= (aW)" 1 =$ 0J B(#)" 1 (20) 
Using Eq. (18), the explicit expression for p^ can be written down: 

mB,e)= p ^^M= P % { B,B c2{ m (2D 

This is ct second important results of this work: similarly to tj c , p v '' retains the same anisotropy of the flux flow and 
pinning tensors, given by the mass anisotropy tensor M alone. Moreover, it satisfies the angular scaling law (with 

a BGL scaling factor 23 s p = 1, as 2,23 pg and pc)- It is important to remember that pf 1 of Eq. (21) is an intrinsic 
quantity, and as such not necessarily it is directly measurable. The actually measurable tensor, according to the 
results developed for the flux flow regime in Ref. 2, is: 



p v {B,9,4>)=p%{B/B c2 {9)) 

with M as: 2 



Ml 



(22) 



M(9,cf>) = -B x (9,cj>)M- 1 B x (9,cf>) (23) 

and where: 

e 2 (9) = cos 2 9 + 7~ 2 sin 2 9 (24) 

is the square of the well-known angular-dependent anisotropy parameter 24 ' 26 which defines, among the others, the 
angle dependence of the critical field B c2 (9) = B c2 (0)/e(9). From Eq. (22) it can be noted that, contrary to the 

=(i) = 

intrinsic resistivity tensor p v , the measurable tensor p v does not obey the angular scaling law, since it incorporates 
additional angular dependencies through M. and e 2 . 



6 



Now I include the effects of flux creep. The pinning energy U depends only on the magnetic field magnitude and 
direction and not on the direction of vortex motion. Moreover, it obeys the usual scaling law (with a constant scaling 
factor 23 7 _1 ). Therefore the thermal depinning time Ttn (Eq. (2)) is a scalar u) v = t~ 1 . The pinning constant intrinsic 
tensor (Eq. (17b)) can thus be corrected to include creep as: 

Consequently the scalar vortex motion resistivity within flux creep, Eq. (5), can be generalized to the anisotropic 
case as follows, yielding an expression analogous to Eq. (21): 

P^{B,9) = pWll^f M = p ^B/B c2 (0))M (26) 

Wo 

where the characteristic frequency a>o and the creep factor e' remain scalar values as in the isotropic case, reported 
in Eqs. (6) and (7). This property will prove important in the interpretation of the experiments. 

Finally, the remaining quantity to be derived is the complex mobility tensor p, V} which is of particular interest for 
its application within the full electrodynamic model to be treated in Section IV. Since the quantity which appears in 
the electrodynamics expressions is actually (B x jl v B x ), 2 it is sufficient to consider the following implicit expression of 
fl v in terms of p = v : 2 

p v = -B x {p\ v $ B)B x (27) 
Incidentally, it is interesting to observe that by comparing Eq. (27) to Eq. (22), one has: 

B x (n v $oB)B x = B x (^^M-^j B x (28) 

The above identity does not allow to equate the tensor between round parentheses because the matrix B x is not 

invcrtible. Actually it can be easily verified that (fi v &oB) ^ (Pff\i/ e2 M" 1 ). Nevertheless, Eq. (28) shows that 
they are physically equivalent when they are involved in the electrodynamics expressions, which always include 
the "sandwich" product with B x . Indeed, this fact can be visually shown by expressing Eq. (28) in a rotated 
frame of reference 2 where the z-axis coincides with the magnetic field direction B. In this frame, the two matrices 
under consideration become identical apart from the elements belonging to the third row and to the third column. 
Considering that these elements do not affect the "sandwich" product with the rotated B x , the equivalence from a 
physical electrodynamics point of view of the two matrices is understood. 

Before concluding this Subsection, it is worth to stress that the choice of a relaxation of the pinning constant 7 to 
deal with flux creep is not a limiting factor to the results obtained up to now: in fact, any pure thermal creep process 
(independent on the angle between the field and the current), possibly with a different definition of e' and uj (Ref. 
14), would yield the same results. 



III. APPLICATION TO EXPERIMENTS: THE MEASURED COMPLEX VORTEX RESISTIVITY IN 

COMMON SETUPS 

In this Section I consider explicitly some typical experimental configurations and I derive specific expressions relating 
measured and intrinsic properties. 



1. Straight planar currents 
The resistivity measured along the current direction J can be computed as: 2 



P (J) 



(pj) ■ J (29) 
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By applying a straight a.c. current J || x and using the resistivity tensor given by Eq. (22), one obtains for the 
measured vortex resistivity: 

pl x \B,6,4>) = P %{B/B c2 (9))f L ^(9^) (30a) 

„ , n ,n 7^ 2 sin 2 sin 2 </> + cos 2 9 s 

# = -2 ■ 2n7 2a ( 30b ) 

7 ^ sin & + cos^ 9 

I recall that <j> is the angle between the projection of the B field on the x-y plane and the x axis (Fig. 1). It can be 
seen that the effective, measured vortex motion resistivity consists in the product of two terms: the first one relates 

(i) 

to the intrinsic resistivity p v n only, which in particular obeys the angular scaling law; the second one, denoted in the 
equation as /l<p(9, </>), contains an "extrinsic" angular dependence which arises from the Faraday and Lorentz actions. 
Consequently, as already anticipated commenting the whole tensor p v , the experimentally measured quantity does 
not obey the scaling law. Only in the case <f) = 7T /2 (B _L J) one has direct, experimental access to the intrinsic vortex 
resistivity, which obeys the angular scaling law. 



2. Rotational symmetric planar currents 



A rotational symmetric planar geometry is often used for measurements of the vortex-state microwave response. 
Examples are cylindrical resonators, 17 ' 27 in which (using cavity perturbation techniques 28 ) the superconducting sample 
is located on the circular bases, and the so-called Corbino disk setup, 9 ' 29 in which the superconducting sample short- 
circuits an open-ended coaxial cable. In both cases, rotational symmetric currents (circular and radial for the resonator 
and the Corbino disk, respectively) are induced along the superconductor a-b plane (see Fig. 2), in the frequent case 
in which the superconductor c-axis is perpendicular to its surface. The effective measured 30 ' 31 resistivity (see the 
following section for the relationship between measured surface impedance and the resistivity) comes out as an 
angular average 32 over the current pattern: 



pi° ] (B,9) = -L jf (MBAfiAaj) ■ J(a)da = 
J (M(9,4>)J(a)Y J(a)da 



1 Pv,u(B,9) 



(31) 



2tt e 2 (9) 

where a is the angle between the (local) current density J and the x axis (see Fig. 2). Exploiting the uniaxial 




a* ( 

1 

-superconductor 



z 
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FIG. 2. Circular (left panel) and radial (right panel) symmetric current patterns, a is the angle between the (local) current 
density J and the x axis. 

anisotropy together with the circular symmetry of the current pattern, the computation can be equivalently and more 
simply done by averaging Eq. (30) over all possible values of 0: 



The result is: 



pi°\B,B)= P y il (B/B c2 (9))f L (9) 
h(0) 



^7 2 sin 2 i 



cos 



7 2 sin 9 + cos 2 9 



(32a) 
(32b) 
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It can be noted that there is no more the angular dependence on <p, and that the function is a particularization 

of the above introduced f^. Additionally, the same comments proposed for Eq. (30) hold also here. Moreover, 

one can note that, in the present case, if p„\i oc (B / B C 2{6)Y , Eq. (32) still yields a scaling law. 31 ' 33 However, in 
the interpretation of the experiments one has to be careful and not confuse this artificial scaling function with the 
theoretical scaling expression. 



3. Angle- dependent effective quantities 



It is interesting to note that a typical GR model analysis in a isotropic superconductor extracts the vortex parameters 
Pff, r and k p from the complex measured p v of Eq. (8) as follows: 



Pff = M(Pv) 

9f(p») 



$Hpv) 



r = 



kri 



^Pv) 

r 



Pff 



W(p v ) + &(p v ) 



On the other hand, when performing measurements on an anisotropic superconductor probed with the circular current 
pattern leading to Eq. (32) (apart from the additional </>- dependence, the same holds obviously for the straight current 
setup related to Eq. (30)), this computation would yield the following effective quantities: 



Pff,tf(B,6)=p%] 11 (B/B c2 (0))f L (e) 
r eff {B,0)=r{B/B c2 {6)) 

k<jUB/B c2 (0)) 



k p , eff (B,0) = 



h(9) 



(33a) 
(33b) 

(33c) 



It can be seen that the parameter r — lu p /lu is directly obtained from the measured quantities: this is an interesting 
result, which allows a direct evaluation of the material anisotropy of the system without the need to deal with Lorentz- 
dependent contribution On the other hand, both pg^ff and fc p , e jj show an additional angular dependence 

through Therefore, in the analysis of angular data care must be devoted in correctly extracting the intrinsic 

quantities instead of the effective ones: this requires to evaluate in some way the function, which in turn requires 

the knowledge of the anisotropy factor 7. 

Further comments can be done considering a scaling analysis performed starting from the effective quantities of Eqs. 
(33). Once the intrinsic quantities Pjf u and r are extracted, they can be checked against the scaling prescription. 



If Pff\i is found to satisfy the scaling, and in the same time r or, equivalently Pq'h, is not, this result would 
unambiguously indicate the presence of directional pinning contributions such as extended defects. 



IV. HIGH FREQUENCY ELECTRODYNAMIC RESPONSE IN ANISOTROPIC SUPERCONDUCTORS 

Having discussed the angular dependent vortex motion resistivity, I now include the effect of supcrfluid and quasi- 
particles in the overall electromagnetic response. In the high frequency regime, an electromagnetic wave impinging 
on a superconductor, driven in the mixed state by a static magnetic field B, determines an electromagnetic response 
which is not only dictated by the motion of the vortices but arises from the coupling between the latter and the high 
frequency currents, which include both the normal and the superconducting components. 8,34 

Moreover, the actual physical quantity that can be directly measured is not the mixed state complex resistivity, but 
instead the superconductor surface impedance 16-18 ' 35 ' 36 Z, defined as the ratio between the components of the electric 
and magnetic fields tangential to the separation surface between the superconductor and the outer medium (typically 
vacuum or air). While the complex resistivity, in the local response regime here considered, can be computed as the 
local relation between the electric field and the current density, the evaluation of the surface impedance Z typically 
requires the determination of the electromagnetic field distribution within the whole superconductor. In tilted static 
magnetic fields the computation is far from trivial. Coffey and Clem (CC) 5 ' 6 ' 8 ' 20,21 used an isotropic description of 
vortex dynamics, with scalar vortex mobility, to calculate the response of isotropic superconductors in parallel 8 and 
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oblique static magnetic fields, 20 ' 21 as well as of anisotropic superconductors in arbitrary oriented magnetic fields. 5 ' 6 In 
the following, (i) I calculate the anisotropic, angular dependent response to electromagnetic fields in the mixed state 
by including the full tensor model for the vortex dynamics developed in the previous Section (Section IV A) and (ii) I 
provide an example of application by computing the surface impedance tensor (whose definition is recalled in Section 
IV B) of an anisotropic superconductor in the commonly found thin film geometry 18 ' 35 (Section IV C). 



A. The fully coupled electrodynamics anisotropic model 

The coupling of vortex motion, superfluid and quasiparticles has been addressed, among the others, 8,34 by Coffey 
and Clem. The core of the CC model is a partial differential equation whose solution yields the high-frequency 
part b of the magnetic field induction existing within the superconductor. Within the linear response regime, for an 
anisotropic superconductor in the e lut sinusoidal regime, under the effect of a static and uniform magnetic field B 
with rigid vortices, the equation in the tensor approach here used reads: 



-V x (A 2 (V x &)) + A*oV x (x 2 a nf E?) - b = 

■Vx(Bx^((vxj)xs))) (34) 



— l- 



The uniform static field B and the high frequency b yield the total magnetic induction B tot = B + b. 

The tensors A and a n j are the anisotropic London penetration depth and the normal fluid conductivity, respectively. 
By expressing them in terms of the mass anisotropy only, thus neglecting any scattering time anisotropy and taking 
the Cooper pair mass anisotropy equal to the electronic mass anisotropy 37 M, one has: 

A 2 = X 2 n M (35a) 

O n f = Onf,\\M~ X (35b) 

where, with the same notation of Eq. (11), X\ Y and <7 n f t n denote respectively the squared London penetration depth 
and the normal-fluid conductivity in the a-b plane. Then, with Eqs. (28) and (22) for p, v and p v , Eq. (34) becomes: 

0b- A^V x (m(V x = 
4<n Vx (^(( Vxg ) x ^)) ) ( 36 ) 



where 

0={-2iX 2 J5 2 nf ^-l) (37) 

and the quantities <5 n /,n and 5 v> u are respectively the normal fluid and vortex motion penetration depths in the a-b 
plane: 

Slf,u = (38a) 

<n = ^ (38b) 

Equation (36), by including the tensor description of the vortex dynamics p, v and p v developed in the previous Section, 
is a very general formulation of the mixed state electromagnetic response, and it includes the effects of the anisotropy 
and of the angle between B and J (see Eqs. (22) and (27)). 

The high frequency field b can be completely determined through Eq. (36) with the appropriate boundary 
conditions. 6 ' 21 Other quantities of interest are the total (superconducting plus normal fluid) current density, which 
can be derived by using the fourth Maxwell equation V x b = /iqJ (displacement current is neglected as ordinarily 
done in superconductors 6 ), and the high-frequency electric field E, which can be derived from the London equation 
yielding: 

E = p cc J (39) 
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where p cc is the Coffey-Clem-like resistivity tensor: 6 



1 , 2iA?! 
1 + 72 



(40) 



where the vortex resistivity in the form of the tensor p v has been explicitly included. 

A significant limit of the above expression is obtained for temperatures T — > 0: in this case a n f,i\ 0, #n/,n ~~ * 00 > 
so that: 



(*) 

i^oA^M + ^i^B.M-^,) (41) 



which shows that at low temperatures the whole complex resistivity is dominated by vortex motion and supercurrents 
and which is the anisotropic variant of the isotropic Brandt's full model. 7 

Once both b and E are known, the surface impedance tensor Z can be computed following the definition valid for 
anisotropic media recalled in the next Section. 



B. Surface impedance for anisotropic media 

For an anisotropic media the surface impedance Z is a two-dimensional tensor defined through the following 
expression: 38 

E„= z(nxh^ (42) 

where n is the normal surface enclosing the media and pointing outwards, toward the surrounding vacuum or air, h 
is the magnetic field equal to bj po since p r = 1, and the subscript "//' denotes the vector components parallel to the 
surface. Placing the superconductor in the positive z-semi-space, its delimiting surface is at z — and in the above 
equation one has n = —z, E ff = (E x ,E y ) and h„ = (h x ,h y ). Using the fourth Maxwell equation V x h = J along with 
E = pJ, one writes: 



E = pVxh (43) 



Combining Eqs. (42) and (43), one obtains: 



Pll Pl2 
P21 P22 




(44) 



From this equation it is clear that in order to compute the surface impedance tensor both the superconductor resistivity 
tensor and the electromagnetic field distribution in the superconductor are needed. 



C. Application to experiments: surface impedance for superconducting anisotropic thin film in titled field 

In this section the expression for the surface impedance of an anisotropic superconductor in tilted magnetic field 
within the CC model, including the anisotropic vortex dynamics description developed in Section II, will be derived 
for a thin film geometry, which is encountered in a large number of experimental studies. 18,35 

Superconducting thin films grown on dielectric substrates are partially transparent to the electromagnetic radiation, 
so that in microwave measurements they are customarily backed with a metallic plate. The resulting layered structure 
is depicted in Fig. 3. The structure is considered to be indefinitely extended along the x-y directions, thus neglecting 
border effects. 

Considering an impinging TEM wave propagating on the z direction, the high frequency fields are uniform on the 
x-y planes. The surface impedance tensor Z can be computed by applying the model presented in Sec. IV A. The 
full computation is provided in the Appendix A3. 

Here I report the result: within the well-known thin film approximation 39,40 t s -C min(<5„/,A) (being t s the film 
thickness), one has: 

2 _ _L | Pcc, 11 Pcc, 12 
t s \Pcc, 21 Pcc, 22 
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FIG. 3. Thin film, grown on dielectric substrate, with backing metal plate. 



which appears as an intuitive (but not straightforward) generalization of the corresponding isotropic expression Z = 

Pcc/ts- 

For the thin film geometry here considered, it can be noted that the surface impedance tensor Z is directly propor- 
tional to the resistivity (sub)tensor. In particular, this direct proportionality allows for the immediate generalization 
to Z of the results concerning the measurable resistivities discussed in Section III, as well as a direct separation of 
the real and imaginary parts of p in the real and imaginary parts of Z. On the other hand, for thicker films the 
expression becomes quickly more cumbersome: by examining the next higher order terms reported in Eq. (A10), it 
can be seen that the Z elements do not correspond to single elements of p cc . Instead the elements of Z become a 
function of the elements of p cc m through complex coefficients, making more difficult the extraction of the resistivity 
from the measured surface impedance. 



V. SUMMARY 

In this paper a full tensor model for the a.c. vortex motion resistivity, including creep, pinning and flux flow, has 
been presented. The model has been employed to develop a full tensor theory for the electromagnetic response of an 
uniaxially anisotropic type-II superconductors in the vortex state. Arbitrary orientations between the applied field, 
the applied current and the anisotropy axis have been considered. Limiting cases have been commented. Relations 
between measurable quantities and intrinsic material properties have been given, showing that care must be put 
in separating the material angular dependence from the one arising from the geometry of the setup. Finally, the 
expression of the measured surface impedance in the largely used thin film geometry has been computed. 
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Appendix A: Application to anisotropic thin films 

In this Appendix the computation of the surface impedance tensor of the thin film geometry described in Sec. IV C 
is fully developed, along the lines of Refs. 5, 6, and 21, where similar computations for different geometries were 
performed. For the ease of comparison, whenever possible the same notation will be used. Moreover, in the following 
the subscript "11" for the a-b plane quantities An and S Vj u is dropped in order to simplify the notation. 
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1. Geometry of the problem 

The geometry here considered is represented in Fig. 3, with the z axis normal to the superconductor surface z = 
and pointing inward. The impinging TEM electromagnetic field is assumed to determine a field: 

b(z = 0) = b y (Al) 

(with bo <C Bo) on the superconductor surface, thus defining the first boundary condition. 
Within the superconductor of thickness t s , i.e. for < z < t s , the field will be in general: 

b{z) = h{z)x + f 2 {z)y (A2) 

Concerning the second boundary, typically the (effective) surface impedance of the film substrate is much larger than 
the characteristic surface impedance of the superconductor, 39 ' 40 a part from very specific cases in which the substrate 
give rise to more complex phenomena. 40 ' 41 At the interface between the superconductor and the substrate, i.e. at 
z = t s , the reflection coefficient for the electric field is ~ 1. Hence the tangential component of b is ~ at the 
interface. 42 

Hence, the overall boundary conditions are the following: 

/i(0) = 0,/ 2 (0) = 6o (A3a) 
h(t s )=0,f 2 (t s )=0 (A3b) 

2. High frequency magnetic field computation 

Applying Eq. (36) to this configuration and separating the two components of b, one has: 

afi (z) + P fi (z) + c/a (z) = (A4a) 
4fk*)+PM*) + cfi(z)=0 (A4b) 
where /3 is defined in Eq. (37) and the parameters a, d and c are defined as: 

d=X 2 - 1 -6 2 v ( 7 - 2 B 2 y + B 2 z )^ (Aoh) 

yh- 2 B xBy ± 

As it can be seen, the parameters a, c and d depend on the superconductor properties and on the static magnetic 
field orientation, here written in compact way as B = (B x , B y , B z ). Equations (A4) and Eqs. (A5), written for the 
thin film geometry, are similar to those reported in Refs. 5 and 6. In Eqs. (A4) and (A5) the vortex mobility tensor 
has been used. 

The coupled differential Eqs. (A4), with the boundary conditions, Eqs. (A3), can be solved by means of the Laplace 
transformation Fi(s) — Jq 00 f i (z)e~ sz dz and anti-transforming the solutions using the partial fraction expansion. One 
obtains fi(z) and fo(z) as linear combinations of the four exponentials e ±z / A± with coefficients Aj_k (with j = 1..4 
and fc = 1, 2) which arise from the fraction expansion: 

f k (z)=A hk e z / x -+A 2 , k e- z / x - + 

+ A 3 , k e z/X + + A 4 , fc e- Z/A + (A6) 

The expressions of the coefficients Aj :k , not reported here, include the first derivatives /{(0) and /^(O) arising from 
the Laplace transformation. The pair of complex penetration depths A± is given by the solutions s 2 = —\± 2 of the 
quadratic equation (ad — c 2 )s 4 + 0(a + d)s 2 + fi 2 = 0. Their explicit expressions are: 



a ^;A- ' ! /;; • Ilh- (Arm) 



c=--5h~ 2 B x B y - (A5c) 
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In the well known 39 ' 40 thin film approximation, t s <C min(<5„/, A), one easily obtains |i s /A±| <C \\ff3t s /\\ <C 1. Within 
this condition, the exponentials in Eq. (A6) can be Taylor-expanded and truncated to the first terms, allowing also 
to easily compute /{ (0) and f' 2 (0) through the application of the boundary condition (A3b) : 



fi{z) = b 



cm 



2{ad~c 2 ) 



1+* 



z 

apt 2 



z 

7! 



t s \2(ad-c 2 ) 



- 1 



t 2 s 2(ad-c 2 ) 



(A8a) 
(A8b) 



Truncation to the second order term avoids f\{z) — 0, a too crude approximation. 

The explicit values of the derivatives at z = 0, needed for the computation of the surface impedance according to 
the definition (44), are reported here: 



/i(0) = -b 

im = -b 



cf3t s 



2(ad-c 2 ) 
1 a/3t s 
T s 2(ad-c 2 ) 



(A9a) 
(A9b) 



3. Surface impedance computation 

It is now possible to explicitly compute the surface impedance elements, applying the definition Eq. (44) with the 
resistivity tensor given by Eq. (40) and with the results of Eqs. (A8) and (A9): 

ry PCC,U Pt s (ap cc ^l + Cp cc 12) ^ PCC, 11 , A1 „ s 

Zu ~ 2(ad-c*) ' = (A10a) 

ry Pcc,2\ Pts{ap C c,21 + CPcc.22) ^ Pcc,2\ , l1nU 

221 =— 2(ad-c 2 ) ' = — (A10b) 

where in the last approximate equalities the higher order terms of the Taylor expansion have been finally neglected. 

Now, by solving again the differential equations with b(z — 0) = box, one obtains the expressions equivalent to 
(A10) for the elements Z i2 and Z22, so that it is finally demonstrated that: 



1 



Pec, 11 Pcc,12 



t s \Pcc,21 Pec, 22 



(All) 
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